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ABSTRACT 

We studied the adsorption of a charged protein onto an oppositely charged membrane, com- 
posed of mobile phospholipids of differing valence, using a statistical-thermodynamical ap- 
proach. A two-block model was employed, one block corresponding to the protein-affected 
region on the membrane, referred to as the adsorption domain, and the other to the unaf- 
fected remainder of the membrane. We calculated the protein-induced lipid rearrangement 
in the adsorption domain as arising from the interplay between the electrostatic interactions 
in the system and the mixing entropy of the lipids. Equating the electrochemical potentials 
of the lipids in the two blocks yields an expression for the relations among the various lipid 
fractions in the adsorption domain, indicating a sensitive dependence of lipid fraction on 
valence. This expression is a result of the two-block picture but does not depend on further 
details of the protein-membrane interaction. We subsequently calculated the lipid fractions 
themselves using the Poisson-Boltzmann theory. We examined the dependence of lipid en- 
richment, i.e., the ratio between the lipid fractions inside and outside the adsorption domain, 
on various parameters such as ionic strength and lipid valence. Maximum enrichment was 
found for lipid valence of about (—3) to (—4) in physiological conditions. Our results are in 
qualitative agreement with recent experimental studies on the interactions between peptides 
having a domain of basic residues and membranes containing a small fraction of the polyva- 
lent phosphatidylinositol 4,5-bisphosphate (PIP 2 ). This study provides theoretical support 
for the suggestion that proteins adsorbed onto membranes through a cluster of basic residues 
may sequester PIP2 and other polyvalent lipids. 
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INTRODUCTION 



Some membrane-associated proteins are known to bind to membranes non-specifically through 
electrostatic interactions (Murray et al., 1997; Resh et al., 1999; McLaughlin et al., 2002; 
Murray et al., 2002). These interactions result from the attraction between a cluster of 
charged residues in the protein and the oppositely charged membrane lipids. 

As the charged protein approaches the membrane, it changes the local membrane com- 
position in its vicinity. We refer to this protein-affected region on the membrane as the 
adsorption domain. The lateral fluidity of the membrane allows oppositely charged lipids 
to migrate toward the adsorption domain in order to minimize the interaction free energy. 
Evidence for such redistribution was reported in experimental (Heimburg et al., 1999; Rauch 
et al., 2002; Wang et al., 2002) and theoretical (Harries et al, 1998; May et al., 2000; Fleck 
et al., 2002; May et al., 2002) studies. 

Lipid redistribution was observed also in monolayers (Lee et al., 1994; Lee and McConnell, 
1995) and bilayers (Groves et al., 1997; Groves et al., 1998) that were subjected to external 
electric fields. 

Local changes in lipid concentration may have biological significance. For example, PIP2, 
a polyvalent phospholipid with valence in the range —3 to —5 (Toner et al., 1988; McLaugh- 
lin et al., 2002; Rauch et al., 2002; Wang et al., 2002) participates in signal transduction 
(Czech, 2000; Payrastre et al., 2001; Simonsen et al., 2001). Its average fraction in plasma 
membranes is very low, about 1%, and it is known to be concentrated in specific regions 
of the membrane (Liu et al., 1998; Stauffer et al., 1998). The PIP 2 lipid serves as a sub- 
strate for Phospholipase C (PLC), which cleaves it to two secondary messengers (Katan and 
Williams, 1997). Another component, the Myristoylated Alanine-rich C Kinase Substrate 
protein (MARCKS), containing an amino acid segment of 13 basic and no acidic residues 
(Blackshear, 1993; McLaughlin and Aderem 1995), is believed to form a PIP2 'reservoir' in 
its adsorption domain. As long as PIP2 is concentrated in the MARCKS adsorption domain, 
PLC is inhibited and cannot catalyze the PIP2 hydrolysis (Wang et al., 2001). It is assumed 
that upon demand, by phosphorylating the MARCKS effector segment, these lipids are freed 
for signaling (McLaughlin et al., 2002). Thus the capability of MARCKS to sequester PIP 2 
potentially affects intracellular signaling. 

Lipids of various valences are attracted to the adsorption domain to different extents. 
This electrostatically- induced enrichment is partially balanced by entropy effects that favor 
homogeneous lipid distribution. Theoretical studies of bilayers composed of neutral and 
monovalent lipids, where the lipid mobility was taken into account, showed that the forma- 
tion of a charged lipid domain due to the adsorbed protein is energetically favorable and 
outweighs entropy effects (e.g., Heimburg et al., 1999; May et al., 2000). Recently, Fleck 
et al. (2002) presented a detailed formulation for the interactions between charged objects 
and a membrane composed of lipids of various valences. These studies demonstrate the 
importance of lipid redistribution in the thermodynamics of protein-membrane adsorption. 

Based on the findings of these theoretical studies, we focus in this work on a simplified 
model for the redistribution of different-valence lipids in the adsorption domain of a charged 
protein. We start by deriving an expression that relates the fractions of the various lipids in 
the adsorption domain to their values in the unperturbed membrane. The general expression 
is restricted, however, to relations between concentrations of different lipids. To obtain the 
actual concentration values for each lipid type, we use the Poisson-Boltzmann (PB) theory 
(e.g., Andelman, 1995). The model is then applied to peptide segments, such as a polylysine 
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Figure 1: Schematic view of the protein-membrane system. The protein region interacting with the 
membrane is modeled as a planar surface of charge density ct p , hovering parallel to the membrane at a 
distance h. The unperturbed membrane is composed of various lipids, each having valence and mole 
fraction fa . The interaction region on the membrane, the adsorption domain, has area A and lipid fractions 
ipi. The whole system is embedded in an ionic solution characterized by a dielectric constant e, Debye 
screening length k _1 and temperature T. 

chain, interacting with a membrane composed of uncharged, monovalent, and trivalent lipids, 
corresponding to zwiterionic PC, PS, and PIP2, respectively. Finally, the model is evaluated 
and the biological implications of its results are discussed. 



A schematic view of the protein-membrane system under study is provided in Fig. d In 
the model, the membrane is considered to be an infinite surface composed of k phospholipid 
species. Each phospholipid type i is ascribed a fraction (in the unperturbed membrane) fa 
and a fixed valence Z\. (We do not consider pH-dependent dissociation of charged lipids, 
which was found to have a minor effect compared to lipid mobility; Fleck et al., 2002.) The 
indices % = and % — 1 are assigned to neutral lipids (z = 0) and monovalent anionic ones 
{z\ = —1), respectively, which are always present in biological membranes. For simplicity, 
all phospholipids are ascribed the same headgroup area a. (In a more detailed model, one 
can treat different headgroup areas; Andelman et al., 1994.) 

We assume that the adsorption domain is a "patch" of finite area A and uniform charge 
density cr m , whose size is much larger than the Debye screening length of the solution i.e., 
n 2 A ^> 1. The screening provided by the surrounding ionic solution ensures a cutoff for the 
effect of the adsorbed protein on the membrane, justifying the finite-area assumption. The 
uniform charge density, employed merely for simplicity, can be thought of as an effective or 
average domain charge density. We thus neglect effects related to charge discreteness. (This 
assumption becomes invalid in certain circumstances; we shall return to it in the Discussion.) 
The membrane outside the domain serves as a large reservoir, assumed to be unaffected by 
the protein. Thus, in this study (similar to the description of May et al., 2002, who treated 
many proteins) we regard the membrane as composed of two blocks, the protein-affected 
adsorption domain, with lipid fractions ipi, and the unaffected remainder of the membrane, 
with fractions fa. 

Both blocks are assumed to contain enough molecules to be considered, to a good ap- 
proximation, as macroscopic phases. (A treatment of finite-size effects can be found in May 
et al., 2002.) The charge densities a and <r m , in the protein-free and protein-bound regions, 
respectively, are defined as: 
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where e is the elementary charge. 

The free energy per unit area of the bare (protein-free) membrane, F^°\ is 

F(°) = ^5>ln^ + F£)(<x) (2) 

i 

The first term is the mixing entropy contribution, where T is the temperature in energy 
units (taking the Boltzmann constant as unity). 

We neglect short-range, non-electrostatic interactions between lipid molecules (May et 
al., 2002) except for excluded- volume effects. The second term accounts for the electrostatic 
contribution. Similarly, the free energy per unit area of the protein-bound domain, F, is 

F= -J^iWi + ^esfcm) (3) 

i 

where F es accounts for the electrostatic interactions among the phospholipids and between 
them and the protein. Note that F^ and F es are functions of 0, and ip i: respectively, only 
via the charge densities a and a m as defined in Eq. [TJ These free energies can be calculated 
using various theories, e.g., the commonly used PB theory (see Appendix II). However, at 
this stage of our formulation, we need not specify the expressions for F^ and F es at all. 

The adsorption domain and the rest of the membrane are at thermodynamic equilibrium, 
thus the electrochemical potentials Hi of each lipid type are equal in the two regions. In 
addition, the membrane incompressibility adds two constraints: 



The electrochemical potential of phospholipid i in the protein-free membrane is 

'II da 



(0) ^ {0) T] <f>i M ,~ 
Hi = a = 1 In — + Z{e (5) 



and similarly, in the adsorption domain, 

9F , dF es 
Hi = a—- = T In — + Zie- — (6) 
oipi Vo oa m 

In Eqs. E]and|nithe dependencies on 0o an d ipo arise from the incompressibility constraint, 
Eq. H| They can be viewed as partial surface pressures exerted by the uncharged species 
(i = 0) due to excluded- volume effects. (Their surface pressure is equal to — Tln0 o and 
—Tin ipo, respectively, in the bare membrane and in the adsorption domain.) We still have 
not specified explicit expressions for F^ and F cs . As a particular example, one may assume 
a mean electric potential (as in PB), having a value \I/( )(0) at the bare membrane, and 
then dF^ ] /da m = ^°\0). If, in addition, we set 0o — l (low membrane charge), then the 
familiar expression for the electrochemical potential is recovered, /i- ^ = Tin (pi + ^-e^^^O) 
(and similarly for //j). Note, however, that the validity of our formulation is more general 
than this specific example. 



Equating [if 3 — Hi we g e ^ 



dF e ( s 0) dF es \ _ 1 ln ^0 O 

da da m / z { 
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Importantly, the left-hand side of Eq. [7| is independent of i. We can therefore compare 
the right-hand side of the equation for a certain species i with the same expression for the 
monovalent species (i — 1). This gives a set of equations relating the enrichment ratios for 
the various lipids, 



We refer to Eq. |H]as the relative enrichment equation. It is a set of (k — 2) equations for 
every i ^ 0,1. The relative enrichment equation is one of the key results of this work. In 
the limit of low charge density, 00,^0 — 1 ; and assuming a mean electric potential having 
the values ^^(0) and ^(O) at the unperturbed membrane and the adsorption domain, 
respectively, Eq. [H] is directly related to a Boltzmann relation, 



As derived above, however, the applicability of Eq. |H] is more general; it is more accurate 
and is valid for highly charged membranes and beyond the mean-field approximation. 

In particular, Eq. El implies that the protein does not perturb the local concentration 
of the neutral lipid, ipo/^o — 1- Yet, as we shall see below, in the biologically relevant 
case where the membrane contains a large fraction of monovalent lipid, the neutral lipid is 
significantly depleted from the adsorption domain. Hence, Eq. |H1 rather than Eq. will be 
used throughout our analysis. 

In the absence of protein, the phospholipid composition would not change and both sides 
of Eq. IHlare trivially equal to 1. Another consequence of Eq. IHlis that perturbation of one 
lipid fraction necessarily entails perturbation in others. In cases where the charge density in 
the adsorption domain increases (in absolute value), |cr m | > |<r|, neutral phospholipids will 
be depleted from the domain to allow the entrance of charged ones, (po/ipo > !• The relative 
enrichment of lipid i, ipi/4>i, is then at least that of the monovalent one raised to the power 



As an example, let us consider the binding of a basic protein to a membrane containing 
uncharged zwiterionic lipids (e.g., PC), monovalent lipids (e.g., PS) and trivalent ones (e.g., 
PIP2). According to Eq. |H1 the polyvalent enrichment ratio will be stronger than that of 
the monovalent fraction by at least a power of Zijz\ = 3. This is a strong effect. If the 
monovalent concentration increases twofold, the trivalent concentration will increase by a 
factor of 2 3 = 8. Similarly, a slight decrease in the negative charge of the adsorption domain 
leads to a significant decrease in polyvalent fraction. 

This entropy-driven enhancement of high-valence ion concentration is an example of a 
more general phenomenon manifest, e.g., in the favorable accumulation of dissolved poly- 
valent ions near charged surfaces or polyelectrolytes (e.g., Rouzina and Bloomfield, 1996). 
Note again that the relative enrichment equation is independent of the specific expressions 
for Fes\a) and F es (a m ). It does not depend explicitly on details such as the distance between 
the protein and membrane, or the protein charge. 

Equation |H1 provides us only with a relation between the different fractions ipi. To cal- 
culate the actual values of ipi we need to derive explicit expressions for the electrochemical 
potentials /4 and //j. To this end, we must introduce details of the protein. It is treated, 
for simplicity, as a flat surface of uniform charge density o~ p located at a distance h paral- 
lel to the membrane (see Fig. Q]). This schematic description may be relevant to proteins 
which have a flat cluster of basic residues facing the membrane at close proximity, while 
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fa ~ i exp[-ez i (^(O) - ¥°\0))/T] 



(9) 
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the rest of the charged residues are further away, screened by the ionic solution. We regard 
the protein-membrane distance h as an external parameter determined by other interactions 
(e.g., desolvation effects), which are not taken into account in our theory. We further as- 
sume that h is small enough {h <C y/A), such that the induced adsorption domain on the 
membrane and the interacting cluster on the protein can be taken to have roughly the same 
area A. 

We apply the commonly used mean-field Poisson-Boltzmann (PB) theory (e.g., Andel- 
man, 1995) to calculate the electrochemical potentials fif^ and /ij. 

The applicability of this theory to solutions containing polyvalent ions is questionable 
(e.g., Netz, 2001). However, in the systems discussed here the polyvalent ions (phospho- 
lipids) are restricted to the membrane and their mole fraction is much smaller than that of 
the monovalent lipids, in accord with biological conditions. The more mobile salt ions in the 
aqueous solution are assumed to be monovalent. Polyvalent ions thus enter the PB calcula- 
tion merely as a (minor) contribution to the membrane surface charge. As such they should 
not induce strong correlations, and the PB theory should be applicable. (An exception, 
where the polyvalent lipid is the majority charge in the membrane and simple electrostatic 
models indeed seem to fail, will be presented in the Discussion.) 

We derive (if* and [i% using three alternative levels of approximation, all of which are 
discussed in detail in Appendix I. 

In presenting the three methods we wish to demonstrate that even a much simplified 
approach, involving minimum computation, still gives useful results. The nonlinear problem 
(NLPB) resulting from the PB theory can only be treated numerically. Subsequently, we 
examine a further approximation where the PB expressions are linearized (LPB). This ap- 
proximation is valid when the electrostatic potential \I/ is much smaller than T/e everywhere 
(Andelman, 1995). Although this condition is not fulfilled in the relevant biological systems, 
the two derivations give similar results 

for reasons that are discussed in Appendix I. The LPB approximation allows us to derive 
analytical expressions for /ij, yet solutions for the various lipid fractions (i.e., for the equations 
/4 = Hi) still cannot be obtained in closed form. 

The fact that, in most relevant systems, the electrostatic interactions dominate over 
entropy effects (May et al., 2000) led us to examine one last approximation, in which the 
derivation is divided into two stages. In the first, we neglect entropy when equating fif^ = 
Hi. This implies that the electrostatic potential is uniform along the membrane (i.e., the 
membrane behaves as a perfect conductor). The resulting membrane charge density in the 
adsorption domain is 

smh(Kh)a - a p 

= t—— (10) 

cosh{Kh) 

In the second stage, this value of a m is substituted in Eq. [T] Equations [H HI and|H] 

(the latter incorporating the effect of entropy) thus provide a closed set of k polynomial 
equations which can be easily solved for the k lipid fractions ipi. We refer to this scheme as 
the simplified linear Poisson-Boltzmann method (SLPB). 

RESULTS 

To study the effects of protein adsorption on a mixed membrane, we calculated the lipid 
fractions for several representative conditions. Our aim is to examine the redistribution of 
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Figure 2: Protein-membrane distance effect. Enrichment ratios (i/>i/</>i) of uncharged (dotted), monovalent 
(dashed) and trivalent (solid) lipid fractions in the adsorption domain of a protein as a function of the protein- 
membrane distance h, scaled by the Debye length The unperturbed membrane contains 69%/30%/l% 
uncharged/monovalent /trivalent lipid fractions. The charge density of the protein is 1.3-fold that of the 
unperturbed membrane and of opposite sign. The inset shows the actual lipid fractions. Although quite 
similar in charge density to the membrane, the protein causes a marked change in lipid composition. Note 
the relatively large increase in the trivalent lipid fraction. 

different-valence lipids in the adsorption domain of a membrane-adsorbed protein as a func- 
tion of several parameters: protein-membrane distance, protein charge, and the valence of 
the most charged lipid species. Unless otherwise stated, we use physiological values for the 
Debye length = 10 A), temperature (300 K), dielectric constant of water (e = 80), and 
lipid headgroup area (a = 70 A 2 ). Throughout the text, the notation 69%/30%/l% un- 
charged/monovalent /polyvalent is used to describe the unperturbed membrane composition 
(4>i values). For consistency, the results presented in this section were all obtained using the 
more elaborate NLPB method. 

Enrichment as a function of protein— membrane distance 

The enrichment ratio, ipi/4>i, was calculated for the association of a charged protein with a 
membrane composed of 69% neutral (zwiterionic), 30% monovalent and 1% trivalent lipids 
(in the absence of a protein). This implies an average charge density a = — 0.33e per lipid 
headgroup area. We present the enrichment arising from a protein that is slightly more 
charged than the membrane, cr p = — 1.3a. (These values are typical to protein-membrane 
interactions, as will be demonstrated later.) Figure El shows the enrichment ratios as a 
function of protein-membrane distance. 

Far from the membrane (h ^> k -1 ) the protein charge is screened and its effect on the 
membrane is weak. As the charged protein approaches the membrane, oppositely charged 
lipids move into the adsorption domain while the neutral lipids are depleted from it. The 
choice of similar charge densities (in absolute value) for the protein and membrane leads 
to minor changes in the fractions of the abundant (neutral and monovalent) lipids, even 
at distances smaller than k~ x . Notably, the fraction of trivalent lipids changes by a much 
larger factor of about 3. This result is a consequence of the exponential dependence of the 
enrichment on lipid valence as seen in Eq. |HJ 

The nonmonotonic behavior at small distances, shown in Fig. is a delicate point that 
deserves further discussion. If the membrane charge density had a fixed value a m ^ — o~ p , 
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then, at a sufficiently short distance, the mutual attraction between the surfaces would turn 
into repulsion (Parsegian and Gingell, 1972). This is caused by the increased concentration 
of the salt ions, which are bound to remain in the confined volume between the protein 
and membrane in order to neutralize the system. In our case, however, the system has the 
additional freedom to change a m . As a result, the electrostatic contribution to the free energy 
of the protein-membrane interaction decreases monotonously with decreasing distance, i.e., 
the interaction is purely attractive (cf. Fig. ITTj). 

(The effect of charge mobility on the interaction between two membranes has been studied 
in detail recently; Russ et al., 2003.) As long as the two objects are not too close, it may 
become favorable to overcharge the membrane and gain attraction energy. This is what 
happens in the system of Fig.|2]for nh < 1. For example, at nh = 0.5, we find |cr m /cr p | ~ 1.17. 
In such a case of overcharging, as the distance is further reduced, the osmotic pressure of 
the salt ions at short distances causes the membrane to decrease its charge density in order 
to lower the energetic penalty of further compression (Fig. |21 in the range nh < 0.5). At 
contact (nh = 0) we have a m = — cr p , such that the system is neutral without mobile ions. 
Thus, the ability to redistribute the lipids allows the system to avoid high concentration of 
ions in the solution. 

Effect of protein charge 

In Fig. El we present the enrichment ratios (V'i/^i) of the different lipid species as a function 
of the protein charge density for a given protein-membrane distance and membrane compo- 
sition. As expected, when the protein is highly charged, the adsorption domain is strongly 
enriched with oppositely charged lipids. On the other hand, at low protein charge, depletion 
of charged lipids is observed. Remarkably, in both cases, the trivalent species exhibits a 
much stronger effect than the monovalent one. This is again a consequence of the sensitive 
dependence of the enrichment ratio on lipid valence (Eq. |Hj). In between the strongly charged 
and weakly charged limits, there is a value of a p for which the membrane is unperturbed 
(see arrow in Fig. EJ). This point does not correspond to |cr p | = \a\, as might have been 
expected. As discussed in the previous subsection, for a nonzero protein-membrane distance 
the membrane may become overcharged. As a result, this special point where ipi/<pi = 1 is 
obtained for |er p | < For example, for the parameters of Fig. [H] it occurs at \a p /a\ ~ 0.57. 

Effect of lipid valence 

Equation |H1 implies that, the higher the valence \z%\ of the lipid, the stronger its enrichment, 
"02/02; relative to that of the monovalent species, ipi/4>i- This does not imply that ^2/^2 V er 
se (not relative to ipi/4>i) increases monotonously with \z 2 \- In fact, there is a competition 
between two opposing effects. The first, which is entropy-driven, favors charging of the 
adsorption domain by high- valence lipids in order to minimize the perturbation of membrane 
composition. On the other hand, from simple stoichiometry, only a small concentration of 
a high-valence lipid is needed to attain a given membrane charge density. The competition 
should result in a maximum of ip2/4>2 at a certain value of valence Z2 = This is confirmed 
in Fig. HJ where we present ip2/<p2 as a function of Z2 for a wide range of protein charge 
densities. (To obtain smoother curves we calculated ip2/ ( t>2 also for artificial, noninteger 
values of z 2 .) For high protein charge, where electrostatic interactions are strong, the entropy 
effect is negligible and \z 2 \ is small; for unphysically high protein charge \z 2 \ eventually tends 
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Figure 3: Protein charge effect. Enrichment ratios (ipi/cfii) of the uncharged (dotted), monovalent (dashed) 
and trivalent (solid) lipids in the adsorption domain of a protein as a function of its charge density (in units 
of elementary charge per lOOOA 2 ). The inset shows the actual lipid fractions. The unperturbed membrane 
contains 69%/30%/l% uncharged/monovalent/ trivalent lipid fractions, corresponding to a = — 4.7e per 
lOOOA 2 and the protein-membrane distance is nh = 0.3. The enrichment in charged lipids increases with 
protein charge. The arrow indicates the tr p value for which the membrane is unperturbed. 



to \z% | = 1. At low protein charge, entropy dominates and \z%\ increases. 

As seen in Fig. EJ \z%\ does not drastically change with protein charge. For reasonable, 
physiological charge densities (cr p of up to ~ 10e/1000A 2 ) we find \z%\ ~ 3 to 4. Thus, lipid 
enrichment due to the adsorption of an oppositely charged protein will be most effective for 
a certain lipid valence which, within our idealized model assumptions, seems to be around 
(-3) to (-4). 

APPLICATION TO PEPTIDE-MEMBRANE INTERACTIONS 

In this section we compare the qualitative results of our model with recent experimental 
studies of the lateral sequestration of the polyvalent lipid PIP2 by adsorbed basic peptides. It 
should be borne in mind that our simple model can only provide an approximate description 
of such systems. Treating the interaction region of the peptide as a large flat surface is a 
particularly severe simplification. We return to this and other weaknesses of the model in 
the Discussion below. 

We focus on two peptides for which there are available experimental data: FA-MARCKS(151- 
175) and a polylysine chain of 13 residues, (Lys)i3. The former corresponds to the basic ef- 
fector segment of the MARCKS protein, where five alanine residues were substituted for the 
original phenylalanine ones (Gambhir et al., manuscript submitted for publication; Wang 
et al., manuscript submitted for publication). We avoid dealing with the MARCKS pep- 
tide, since experiments indicate that its hydrophobic phenylalanine residues pull the peptide 
into the membrane such that its backbone penetrates the membrane (Qin and Cafiso, 1996; 
Zhang et al., 2003). This kind of interaction is expected to be sensitive to specific molecular 
details, which are not encompassed by our model. 

In order to apply the model, we need an estimate for the area of the peptide that interacts 
with the membrane in order to determine the peptide charge density. In addition, this area is 
assumed to be equal to that of the adsorption domain, A (Fig. 1). We built extended peptides 
(MOE software, 2002), similar to the ones used by Wang et al. (manuscript submitted for 
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Figure 4: Valence effect. The enrichment ratio of polyvalent lipids (^2/^2) in the adsorption domain of 
a protein as a function of the lipid valence z-i. The different curves correspond to different protein charges, 
the number on each curve indicating the number of elementary charges per lOOOA 2 . Maximum enrichment 
is marked by a box. The unperturbed membrane composition and protein-membrane distance are the same 
as in Fig. 

publication). We then denned the effective peptide area as the area of its projected backbone 
plus an envelope of width k~ x around it. (This somewhat arbitrary definition will be further 
examined at the end of this section.) Dividing the number of charged residues in the peptide 
by this area, we got the estimates cx p ~ 13e/2120A 2 and 13e/1060A 2 for FA-MARCKS(151- 
175) and (Lys)i 3 , respectively. Note that for these calculated areas and typical peptide- 
membrane distances h of a few angstroms, the basic assumption of the model, h <C VA, is 
well satisfied. 

We used these estimated a p values to produce Fig. plotting the trivalent lipid fraction 
in the adsorption domain of each peptide as a function of its distance from the membrane. 
In accordance with experiments (Gambhir et al., manuscript submitted for publication), we 
took the unperturbed membrane composition to be: 0o — 82% (corresponding to the un- 
charged zwiterionic PC lipid), 4>\ — 17% (monovalent PS), and 2 = 1% (trivalent PIP2). 
Figure shows that the PIP2 fraction rises to 18% in the (Lys)i3 adsorption domain, whereas 
only 6% PIP2 is obtained in the case of FA-MARCKS(151-175). This is caused by the 
higher charge density of (Lys)i3, roughly double that of FA-MARCKS(151-175). However, 
FA-MARCKS(151-175) has twice the effective area of (Lys)i3; thus, if we examine the av- 
erage number of PIP2 molecules per peptide adsorption domain, N 2 , the difference is less 
significant— while N 2 ~ 2.7 for (Lys)i 3 at nh < 0.3, for FA-MARCKS(151-175) at the same 
distance N2 ~ 1.8. (The value of nh ~ 0.3 corresponds to h ~ 3A at lOOmM salt, which 
is the approximate peptide-membrane distance; Ben-Tal et al., 1996; Murray et al., 2002. 
Consequences of this small value will be addressed in the Discussion below.) 

Next, we examined the dependence of the trivalent lipid (PIP2) fraction in the adsorp- 
tion domain on 1; the monovalent lipid fraction in the unperturbed membrane. Figure El 
presents the results obtained for FA-MARCKS(151-175) and (Lys)i3 at various distances. 
The smaller the value of <fii, the stronger the enrichment in trivalent lipids. When there is 
a little amount of monovalent lipids, the membrane charge density induced by the peptide- 
membrane interaction is attained primarily by the trivalent species. Therefore, at <pi = 0, 
the enrichment in trivalent lipid is maximum. At that limit, the number of PIP2 molecules 
per adsorption domain can be simply approximated (at short distances) as the number of 
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Figure 5: Comparison between lipid rearrangement induced by FA-MARCKS(151-175) and (Lys)i3. Triva- 
lent lipid fraction in the adsorption domain of the FA-MARCKS(151-175) (dashed) and (Lys)i3 (solid) 
peptides as a function of their distance from the membrane. The unperturbed membrane is composed of 
82%/17%/l% uncharged/monovalent/trivalent lipid fractions. The enrichment caused by (Lys)i3 is much 
stronger than that achieved by FA-MARCKS(151-175) due to its higher (roughly double) charge density. 

charges on the peptide divided by the lipid valence, N 2 ~ \Aa p /(z 2 e)\. (In practice, how- 
ever, this example of the polyvalent lipid being the majority charge is problematic, as will 
be presented in the Discussion.) Figure |H1 shows that, under physiological conditions of 30% 
monovalent lipid fraction and only 1% PIP2, both peptides sequester PIP2 — roughly one 
molecule per FA-MARCKS(151-175) peptide and around two molecules per (Lys)i3 pep- 
tide. This result is in qualitative agreement with experiments (Gambhir et al., manuscript 
submitted for publication), where (Lys)i3 was found to attract PIP2 more strongly than 



To emphasize the strong sequestration of the trivalent lipid, consider the case of a 
73%/17%/0 PC/PS/PIP2 membrane interacting with (Lys)i3. At a short distance the pep- 
tide charge will be neutralized by about 13 PS molecules. For a 72%/17%/l% membrane, 
as shown by the dashed curve in Fig. |U1 there are roughly 2-3 sequestered PIP2 lipids per 
(Lys)i3, providing 6-9 elementary charges. Thus, introducing one percent PIP2 into a mem- 
brane of 17% PS leads to replacement of about one half of the PS lipids in the adsorption 
domain by PIP2. Indication of such an exchange of PS for PIP2 upon addition of a small 
amount of PIP2 was found in a recent experiment (S. McLaughlin, personal communication). 

Interestingly, the concentration, on average, of about two PIP2 molecules per (Lys)i 3 
would not be possible if PIP2 were of much different valence. Figure [7| shows the average 
number of PIP2 molecules per adsorption domain of both peptides as a function of the 
PIP2 hypothetical valence. Similar to the results shown in Fig. 01 we find a nonmonotonic 
behavior as a function of valence with a maximum at \z 2 \ ~ 3 to 4. It is stressed again 
that, in view of our simplified model, one should pay more attention to the existence of a 
competition mechanism, leading to an optimum valence, than to the exact value obtained 
for that valence. 

In Fig. |S1 we show the dependence of the trivalent lipid fraction in the adsorption do- 
main on its value in the unperturbed membrane. As expected, ip 2 increases with (p 2 . This 
calculation shows that one needs 4> 2 > 1% in order to get an average stoichiometry of 1:1 
between FA-MARCKS(151-175) and PIP2- As we have seen above, (Lys)i3 sequesters PIP2 
more effectively. Hence, as demonstrated in Fig. |SJ a <p 2 value of only 0.1% is sufficient 



FA-MARCKS(151-175). 
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Figure 6: Effect of monovalent lipid concentration. Trivalent lipid fraction (left ordinate) and number 
(right ordinate) in the adsorption domain of (A) FA-MARCKS(151-175) and (B) (Lys) 13 peptides as a 
function of the peptide-membrane distance. Different curves correspond to membranes of different lipid 
compositions: 69%/30%/l% (solid), 72%/17%/l% (dashed), and 99%/0%/l% (dotted). The enrichment in 
the polyvalent lipid increases with decreasing monovalent-lipid fraction. 



to obtain a 1:1 (Lys)i3:PIP2 ratio. That is, (Lys) 13 can sequester an appreciable amount 
of PIP2 even when the membrane contains 300- fold more PS lipids than PIP2. This value 
is in very good agreement, perhaps fortuitously, with the measurements of Gambhir et al. 
(manuscript submitted for publication). 

Finally, we examined the dependence of the PIP2 enrichment on the ionic strength, i.e., 
the concentration of mobile salt ions in the solution, uq. This parameter enters into the 
model through the Debye screening length k; -1 , which both scales the distance h and affects 
the amplitude of the electrostatic potential (see, e.g., Eq. |2T]in Appendix I). Figure M shows 
that ip2 decreases with ionic strength. This is a result of the increased screening of the 
electrostatic attraction between the membrane and protein. The changes are not dramatic 
up to quite high no values. The reason is the very close proximity of the two objects (3A) 
for which nh < 1 in the entire no range examined. 



Effect of approximated peptide size 

As discussed above, we define the effective area of the peptide as the area of its projected 
backbone plus an envelope of width k~ x around it. This definition, however, is somewhat 
arbitrary. We therefore examine the effect of relaxing the effective-area definition on the 
values obtained for N 2 - 

Figure E3 shows the average number of trivalent lipids N 2 per adsorption domain of a 
20-amino-acid-long peptide for a range of 1/A values extending to ±50% of our original 
estimate. In this range the alteration in N 2 is limited to ±1 lipid molecules. 



Interaction free energy 

From the Poisson-Boltzmann theory, as applied to our model, we can calculate the contri- 
bution to the free energy of peptide-membrane association coming from electrostatics and 
entropy. The derivation is given in Appendix II. The results for the case of FA-MARCKS(151- 
175) interacting with a membrane of different compositions are presented in Fig. Some 
of these compositions have also been studied by Gambhir et al. (manuscript submitted for 



12 




Figure 7: Effect of hypothetical PIP2 valence on sequestration. Average number of polyvalent lipids per 
adsorption domain of FA-MARCKS(151-175) and (Lys)i3 as a function of PIP2 valence. Concentration is a 
maximum for \z%\ ~ 3-4. Both calculations were performed using a membrane composition of 69%/30%/l% 
uncharged/monovalent /polyvalent lipid fractions in the unperturbed membrane, and a peptide-membrane 
distance of nh = 0.3. 

publication) and Wang et al. (manuscript submitted for publication). It should be recalled 
that our model does not take into account repulsive interactions, such as the Born desolva- 
tion effect (see Wang et al., manuscript submitted for publication). Lipid demixing effects 
are included in the model, so the local membrane charge density may change due to the 
peptide. As a result, the electrostatic free energy in our model decreases monotonously as 
FA-MARCKS(151-175) approaches the membrane. To get the total free energy of associa- 
tion, one ought to add a repulsive interaction at short distances (not included in our model), 
that would yield a free-energy minimum at a distance of a few angstroms (see, e.g., Wang 
et al., manuscript submitted for publication). Thus, the free-energy values presented in Fig. 
rm are probably more negative than the actual binding free energy. More detailed models 
provide values for the total binding free energy that are roughly one half the contribution 
presented here (Wang et al., manuscript submitted for publication). 

As expected, electrostatic attraction between the peptide and the membrane is propor- 
tional to the charge density of the membrane (Fig. Illj). The interaction free-energy is not 
sensitive to the specific lipid composition of the membrane, but rather to its average charge 
density. This result stems from the minor contribution of mixing entropy to the free energy 
(fourth term in Eq. I2*3*j) and the approximate behavior of the membrane in our model as 
a surface of constant electric potential, which is determined by the average charge density 
a (May et al., 2000). The calculated contribution to the free energy, therefore, is mainly 
the work required to bring a charged object (protein) into such a potential, regardless of 
the lipid redistribution. This may also be the reason why FDPB calculations (e.g. Ben-Tal 
et al., 1996; Murray et al., 2002), although ignoring lipid redistribution, often give a good 
estimate of the binding free energy. 
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Figure 8: Effect of (j> 2 on enrichment. Trivalent lipid fraction (left ordinate) and average number (right 
ordinate) in the adsorption domain of (A) FA-MARCKS(151-175) and (B) (Lys)i3 peptides as a function 
of the peptide-membrane distance. The two curves correspond to membrane compositions of 69%/30%/l% 
PC/PS/PIP2 (solid) and 6! 



0.2 




0.05- 



°0 200 400 600 800 1000 

n (mM) 

Figure 9: Ionic strength effect. Trivalent lipid fraction in the adsorption domain of FA-MARCKS(151- 
175) (dashed) and poly-lysine (solid) peptides as a function of salt concentration. The two curves for each 
peptide were calculated using membrane compositions of 69%/30%/l% (upper), and 82%/17%/l% (lower). 
A peptide-membrane distance of h = 3A was used. 

DISCUSSION 
Model evaluation 

A simplified theoretical model for the effect of an adsorbed charged protein on the distri- 
bution of phospholipids in the membrane has been presented. One of the key theoretical 
results is the relative enrichment equation, Eq. |H1 relating the concentrations of various lipid 
species in the adsorption domain. It should be emphasized that this result is almost model- 
independent. Within a two-block picture, it should hold for any protein-membrane system 
governed by electrostatics, regardless of particular details of the protein, the strength of 
the electrostatic interactions, and the validity of a mean-field (PB) assumption. The main 
physical effect described by this relation is the sensitivity of lipid enrichment to valence, i.e., 
the increased concentration of polyvalent lipids compared to that of monovalent ones in the 
adsorption domain. 

We have demonstrated how simple Poisson-Boltzmann calculations can be added to Eq. 
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Figure 10: Dependence of polyvalent lipid number on estimated peptide area. Average number of trivalent 
lipids (N2) as a function of the inverse area of the adsorption domain I /A. A is also taken as the effective 
peptide area. The arrow marks the 1 j A value according to the definition used throughout this study. Results 
were obtained for a 10 (dotted), 15 (dashed) and 20 (solid) charged residues in a 20-amino-acid-long peptide 
interacting with a 69%/30%/l% uncharged/monovalent /trivalent membrane at a distance of nh = 0.3. 

IH1 to obtain further details of the membranal interaction region. In most of the biologically 
relevant cases, including the peptides discussed here, the peptide-membrane distance h is of 
the order of a few angstroms only (Ben-Tal et a!., 1996; 1997; Murray et al., 2001; 2002). In 
the context of our model, the elaborate Poisson-Boltzmann theory is unnecessary for such a 
thin water layer between the two charged objects. Nevertheless, it should be recalled that, 
at any rate, our calculations converge to the correct membrane charge density at small h 
due to charge neutrality, cr m = — cr p (cf. Fig. H2J). 

The model is focused on the qualitative behavior of electrostatically dominated protein- 
membrane systems. Our description of the protein-membrane interaction is evidently crude, 
as it omits the molecular details of this complex system; see the atomistic description used in 
the companion report of Wang et al. (manuscript submitted for publication). Thus, effects 
related to charge discreteness and three-dimensional structure are ignored. Such effects, for 
example, may play an important role in the membrane adsorption of the MARCKS protein, 

whose backbone is believed to penetrate the membrane (Zhang et al., 2003) and whose 
charges are not evenly spaced. Furthermore, we considered only electrostatic and entropy 
effects while neglecting other interactions (e.g., desolvation). The model is not intended 
to reproduce such data as the total binding free energy or the optimum protein-membrane 
distance (e.g., Ben-Tal et al., 1996; Arbuzova et al., 2000). Similarly, the values that have 
been obtained for the various lipid concentrations should be regarded as rough estimates. 

Treating a peptide as a large, flat, uniformly charged surface is probably the gravest 
simplification of the model, deserving further discussion. The assumption of flatness may 
be applicable to proteins interacting with the membrane through such effector regions as 
those described in Wang et al. (manuscript submitted for publication) and Gambhir et al. 
(manuscript submitted for publication). The assumption regarding the lateral extent of the 
peptide is reasonable in cases where the peptide-membrane distance h is much smaller than 
both lateral dimensions. For example, in the case of the FA-MARCKS(151-175) discussed 
above, we estimated the peptide surface area facing the membrane as a rectangle of dimen- 
sions 106Ax 20A. The smaller lateral dimension is still significantly larger than the typical 
value of h ~ 3A. 
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Figure 11: Contribution from electrostatics and entropy within the smeared charge model to the free 
energy of pcptide-membrane interaction. Free energy in units of T as a function of the peptide-membrane 
distance. The curves correspond to FA-MARCKS(151-175) and membrane compositions of 69%/30%/l% 
(solid), 77%/18%/5% (dash-dotted), 72%/17%/l% (dashed), and 90%/10%/0% (dotted). The vertical 
dotted line indicates a distance of nh — 0.3. The presented contribution to the free energy decreases 
monotonously with distance. The solid and dash-dotted curves correspond to membranes of very different 
lipid compositions but the same charge density (0.33e/a). 

In addition, the fact that this smaller lateral dimension (~ 20A) is comparable to the 
screening length may lead, in principle, to significant finite-size effects. Yet, as already 
noted above, because of the small value of h in the relevant systems, the charge density in 
the adsorption domain should become insensitive to these details and be determined, to a 
good approximation, merely by charge neutrality. 

The small value of h, on the other hand, raises a difficulty with respect to the smeared- 
charge simplification. Since h is similar to or smaller than the typical distance between 
charged groups on the peptide, neglecting charge discreteness is clearly questionable. Spread- 
ing the charges of the lipid headgroups evenly over the membrane is problematic as well. 

Nonetheless, several recent experiments might help us indicate the limits of validity of 
models based on smeared electrostatics, such as ours. The affinity of FA-MARCKS(151- 
175) to PC:PS membranes was found to depend linearly on PS fraction (Murray et al., 
1999). Adding a small amount of PIP 2 to such membranes, having a significant PS fraction, 
did not change the binding significantly (S. McLaughlin, personal communication). These 
two observations are in line with simple electrostatic considerations. However, the binding 
affinity of the same peptide to a membrane of relatively low charge density, composed of 99:1 
PC:PIP2, is surprisingly large — similar to that of a 5:1 PC:PS membrane — even though 
the membrane charge densities in the two cases differ by a factor of about 5 (Wang et al., 
2002). Clearly, the latter observation cannot be accounted for by smeared electrostatics. We 
therefore believe that our model gives reliable results regarding lipid redistribution only in 
cases where the membrane has a high "background" charge density, i.e., a large fraction of 
monovalent lipid, which is the biologically relevant case. For such membranes, as discussed 
in the previous section, the model predictions are in qualitative agreement with available 
experiments. 

Apart from these assumptions, the model contains another implicit simplification, namely, 
that all the electric field lines are contained within the aqueous spacing between the mem- 
brane and protein. This commonly used assumption is strictly correct in the limit where the 
objects are either infinitely thick or of a vanishing dielectric constant. As will be reported 
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elsewhere, we find that this approximation is, in fact, still good for objects of e = 2 and 
thickness as small as one third the Debye length. 

On the positive side, the simplified model presented here provides new insights into 
lipid redistribution caused by protein adsorption. Although most of the results presented 
in this work were obtained using a numerical solution of the nonlinear Poisson-Boltzmann 
equation, we have demonstrated how one can get the same qualitative results using a much 
simpler calculation, the SLPB method, involving merely a set of polynomial equations (see 
Appendix I). Such a scheme may serve as a better starting point for more detailed numerical 
calculations, e.g., FDPB (Honig et al., 1993). In most of the current FDPB calculations, a 
predefined membrane composition identical to that of the bare membrane is used (Ben-Tal et 
al., 1996; 1997; Murray et al., 1999; Wang et al., manuscript submitted for publication). This 
arbitrary description of the adsorption domain may be improved if one uses a preliminary 
analysis of the type presented here to produce an approximate lipid configuration. 

Future extensions of this work may include phenomena such as elastic deformation of the 
membrane (Dan et al., 1993; May, 2000), adsorption of multiple proteins (May et al., 2000; 
May et al., 2002), nonuniform charge density in the adsorption domain (May et al., 2000), 
and acid dissociation at different pH values (Fleck et al., 2002). 

Biological implications 

ENTH, FYVE, PX and other membrane-association domains use predefined stereochem- 
istry to recognize poly-phosphoinositides (McLaughlin et al., 2002; Lemmon, 2003). These 
domains, which are commonly found in proteins involved in intracellular signaling, bind 
tightly to the poly-phosphoinositides, often via ion pairs, anchoring the protein firmly and 
irreversibly to membrane surfaces. The binding specificity is reflected in the evolutionary 
conservation across the homologous domains comprising the family; usually the amino acid 
residues that mediate the poly-phosphoinositides binding are strictly conserved and can often 
serve as sequence signatures to recognize these domains using sequence analysis tools. 

Here we dealt with a much less specific, and often reversible, mode of membrane recog- 
nition via a cluster of basic residues on the membrane-facing region of the protein. These 
residues interact electrostatically with acidic lipids in the adsorption domain on the bilayer 
surface. Our model showed that in such cases membrane association induces a preference 
for polyvalent lipids such as PIP2 to sequester in the adsorption domain. The number 
of sequestered polyvalent lipids may be regulated by the charge density and size of the 
membrane-interaction region on the protein. This result supports the suggestion, advocated 
in the accompanying papers (Gambhir et al., manuscript submitted for publication; Wang 
et al., manuscript submitted for publication), that membrane-associated proteins such as 
adducin (Matsuoka et al., 2000), DAKAP200 (Rossi et al., 1999), GAP43 (Laux et al., 2000), 
MARCKS (Wang et al., 2002), and MacMARCKS (Blackshear, 1993), which contain a clus- 
ter of basic residues, may create a reservoir of PIP2 molecules in their adsorption domain 
(McLaughlin et al., 2002). It may further imply that the membrane-interaction region of 
such proteins does not have to be strictly conserved evolutionarily; it should only preserve 
a specific charge density, as in C2 domains (Murray and Honig, 2002). This speculation, 
naturally, needs to be checked in future studies. 

We found that, when all other parameters are held fixed, there is an optimum value of 
lipid valence that yields maximum enrichment in the adsorption domain (Figs. E]and[7j). 
This value, resulting from a competition between entropy and stoichiometry, is found to be 
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at reasonable valence values (e.g., z* 2 = —3 to —4 in the examples above). Thus, if in a certain 
biological scenario there is a need to increase the local concentration of a phospholipid by 
electrostatic interactions, a polyvalent lipid of valence larger than 1 but not too large (say, 
~ —3) would be advantageous. The valence of PIP2, considered to be between —3 to —5 
(Toner et al, 1988; McLaughlin et al, 2002; Rauch et al, 2002; Wang et al, 2002) appears 
to be in line with this criterion. 

Assuming a certain valence for PIP2 and a certain distance for peptide association with 
the membrane, our model enables the derivation of a set of approximate rules relating 
the number of basic residues on the peptide to the average number of sequestered PIP2 
molecules per peptide. For example, assuming a trivalent PIP2 and association distance of 
about 3A, under physiological conditions, each segment of seven consecutive lysine residues 
of an adsorbed poly-lysine peptide such as (Lys)i 3 (Fig. |J3) would sequester on the av- 
erage approximately one PIP2 molecule when the membrane composition is 69%/30%/l% 
uncharged/monovalent /trivalent, and roughly one and a half PIP2 molecules when the mem- 
brane composition is 82%/17%/l%. These results are in good agreement with the detailed 
calculations reported in the companion report of Wang et al. (cf. Table 2B in Wang et al., 
manuscript submitted for publication). 
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APPENDIX I: METHODS FOR CALCULATING LIPID CON- 
CENTRATIONS 

In this section we present the method for calculating the values of ipi. We present three 
levels of approximation and compare their results. 

In thermodynamic equilibrium the electrochemical potentials of each lipid species in the 
protein- free membrane and in the protein adsorption domain should be equal, = /ij. 
Within a mean-field approximation, this condition can be written as 

^ (0) (0) +Tln^i = z ie tf(0) + Tin ^ (11) 
<Po Wo 

where \l/( )(0) is the mean electrostatic potential at the bare membrane, and ^(O) its value 
at the adsorption domain. Equation^^is actually a set of (k — 1) equations for every species 
i 7^ 0. The incompressibility condition, YliWi = 1> closes a set of k equations for the k 
unknown ipi. To solve these equations we need the surface potentials \l/( )(0) and ^(0)- We 
derive them using the Poisson-Boltzmann theory. Note that \l/(0) depends on the variables 
ipi that determine the charge density of the adsorption domain, a m . 



Nonlinear Poisson-Boltzmann 

To find the surface potentials we need to solve the Poisson-Boltzmann equation, 

d^y 

— = n 2 sinhy (12) 
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with the appropriate boundary conditions. In Eq. H21 y is the local dimensionless potential 
at a distance z from the membrane, y(z) = e$>(z)/T. 
For the bare membrane the boundary conditions are 



d y(0) 



dz 



z=0 



Airea 
eT 



(13) 



and vanishing of the field at z — * oo. This problem, of a single charged plate in an electrolyte, 
is analytically solvable in closed form (Andelman, 1995). The result is 

(14) 

I. (We have 



l/ (0) (0) 



-4tanh 



7 



where 7 is the positive root of the quadratic equation 7 2 + r yKeT/(7re\a\) 
assumed the membrane to be negatively charged.) 

For the protein-membrane system, the boundary conditions are 
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dz 
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dz 
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(15) 



Integrating the Poisson-Boltzmann equation (Eq. IT2|) once while applying both boundary 
conditions, we get 



dy 

dz 



/■ca/2 coshy + C , C 



( 47rea I i 



coshy(O) — coshy(h) 
We then integrate Eq. HH to get 
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dz 
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y(h) 

j,(o) A/cosh ?/ + C/2 



(16) 
(17) 

(18) 



Equations IT7I and ITH1 are solved numerically for the two unknowns y(0) and y(h), the surface 
potentials of the membrane and protein. Finally, the values of y(0) and y^(0) (Eq. I14|l are 
used in Eq. [TT]to calculate the values of 



Linear Poisson-Boltzmann 

The LPB approximation holds when the electrostatic interactions are weak compared to T, 
?/ <C 1. (In fact, the systems relevant to our study are far from this limit; cf. Fig. El) in 
this limit we can linearize the Poisson-Boltzmann equation (Eq. I12J1 . 



dy 



n 2 y (19) 



dz 2 

The surface potential of the unperturbed membrane then takes the form (Evans and Wen- 
nerstrom, 1994), 

y (0 '(0) = (20) 

enl 

For the protein-membrane system we integrate Eq. EI]twice while applying the boundary 
conditions (Eq. ITH|) to obtain (Parsegian and Gingell, 1972) 

Aire 

yW = e.Tsinh(.h) K + a - C ° Sh{Kh)] (21) 
The potential values of Eqs. EU] an d |^ are then used in Eq. El to calculate the ipi values. 



19 



0.15- 




Figure 12: Comparison of the different computational methods. Trivalent lipid fraction in the adsorption 
domain of a highly charged protein (a p = 13e/1000A) as a function of the protein-membrane distance. 
Results were obtained using three methods: nonlinear Poisson-Boltzmann (solid curve), linear Poisson- 
Boltzmann (dashed curve), and the simplified linear Poisson-Boltzmann (dotted curve). At distances nh < 
0.3, where the membrane approaches charge matching, the three methods differ by less than 10%. The 
unperturbed membrane composition is the same as in Fig. |3J 

Simplified linear Poisson-Boltzmann 

In this further approximation the entropy is neglected in /4 and //j, and we are left with a 
uniform surface potential y(0) = y(°>(0). 

Note that the entropy contribution is not necessarily small compared to T. The strong 
enrichment in trivalent lipid, as demonstrated in Figs. El and entails an entropy penalty 
of a few T. Yet, it is still much smaller than the electrostatic contribution, which amounts 
to tens of T (cf. Fig.EJ. 

Within the LPB approximation we can equate expressions l2*Ul and 12 II to obtain the charge 
density in the adsorption domain in closed form as given in Eq. [HI (May et al. 2000). 
Substituting <r m in Eq. [T] yields a set of k equations, Eqs. HUH and|Hl which are easily solved 
for the k lipid fractions t/>$. 

Comparison of the methods 

Here we compare the NLPB, LPB and SLPB calculation methods for a specific example. 
In the limit of weak electrostatic interactions compared to the thermal energy T the LPB 
method should coincide with the NLPB one. To highlight the difference between the methods 
we therefore chose as an example a highly charged protein (a p = 13e/1000A 2 ). The trivalent 
lipid fractions in the adsorption domain, as calculated using the three methods, are plotted 
in Fig. For both long and short distances all three curves match. In the long distance 
limit (k/i ^> 1), the protein-membrane interaction is weak and the linearization of the 
Poisson-Boltzmann equation is valid. At very short distances the membrane, as described 
by all three methods, is forced by charge neutrality to match the charge density of the 
protein (in opposite sign). In the intermediate range the methods differ, yet because of the 
two constraints at k/i < 1 and nh ^> 1 the differences are mild. It is noteworthy that 
at distances nh < 0.3 (typical to protein-membrane adsorption), the calculated ipi values 
obtained using the three methods differ by less than 10%, and the difference in magnitude 
between the charge densities of the protein and membrane are also less than 10%. 
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APPENDIX II: FREE ENERGY CALCULATION 



In this appendix we calculate the free energy of protein-membrane interaction as arising 
from our model. It should be recalled that what is presented here is not the total binding 
free energy but only the contributions due to simple electrostatics and entropy. In particular, 
since charges are allowed to redistribute in the membrane, the interaction is purely attractive 
and the free energy decreases monotonously with distance, reaching its minimum at h — 0. 
This should be contrasted with the FDPB calculations (e.g., Ben-Tal et al., 1996; Wang et 
al., manuscript submitted for publication) which take desolvation into account and yield a 
free-energy minimum at a finite distance. 

The general form of the free energy per unit area within a mean-field approximation is 



F(h) 



e 



8tt V dz 



2 



+ e(n+ — n)^ 



dz 



+ [T (n+ In n + — n + + n_ In n_ — n_ — 2n In n + 2n ) 
Jo 

—fi + (n + — tiq) — //_(n_ — n )}dz 

T . 1 . 

+ - V ipi In - - V inif>i + a P V(h) + <x m #(0) (22) 

The first term in Eq. |221is the energy associated with the mean electric field. The second term 
corresponds to the interaction of the mobile ions with the field, where n + and n_ are the local 
concentrations of the monovalent positive and negative ions. The second integral accounts 
for the ideal entropy of mixing of the mobile ions and their contact with ion reservoirs having 
chemical potentials /i+ and The last line in expression 1221 corresponds to the surface 
energy of the membrane and protein, including mixing entropy of the lipids and electrostatic 
interactions. 

Setting the variations of F with respect to n + ,n_, and \I/ to zero, one properly recovers the 
Boltzmann relations for the mobile ions, the Poisson equation, and, hence, also the Poisson- 
Boltzmann equation. Using these results along with Eq. and changing the integration 
variable from z to y, we obtain the following simplified expression for the free energy: 



2Tn f^ h H -2coshy-C/2 ^ Ta p (h ,Ta m T^ I ^ 

F{h) = / = — dy + — £y(h) + y (0) + ->VilnV><--> Hi^i 

« J y (o) V 2 cosh y + C e e a ^ a 

" (23) 

where the constant C was defined in Eq. [H3 We can now use the values obtained for y(0), 
y(h), fa, and ipi, as described in appendix I, to calculate the free-energy of interaction F(h). 
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